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Abstract 

We explore the Higgs sector of the Minimal Left-Right (LR) Model based on the gauge group 
SU(2)i x SU(2)r x U(1)b~l with explicit CP violation in the Higgs potential. Since flavour¬ 
changing neutral current experiments and the small scale of neutrino masses both place stringent 
constraints on the Higgs potential, we seek to determine whether minima of the Higgs potential 
exist that are consistent with current experimental bounds. We focus on the case in which the right- 
handed symmetry-breaking scale is only “moderately” large, of order 15-50 TeV. Unlike the case 
in which the Higgs potential is CP-invariant, the CP noninvariant case does yield viable scenarios, 
although these require a small amount of fine-tuning. We consider a LR model supplemented by 
an additional 17(1) horizontal symmetry, which results in a Higgs sector consistent with current 
experimental constraints and a realistic spectrum of neutrino masses. 
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I. INTRODUCTION 


Left-Right (LR) models have long 
dard Model (SM) of particle physics 


provided attractiv 

QBIiiiQ- 


provided attractive extensions to the well-tested Stan- 
A typical LR model is constructed 


by enlarging the SM gauge group SU (2) L xU (l)y to the group SU (2) L x SU (2) R xU(l) b-l 
and thus contains an extra set of W and Z gauge bosons compared to the Standard Model, 
as well as an extended Higgs structure. Minimal versions of the left-right model usually 
contain one bidoublet Higgs boson, 0, as well as left- and right-handed triplet Higgs bosons, 
A l,Ri although other constructions are also possible (see for example Ref. jg]). 

There are several reasons for the appeal of LR models. Among these is the fact that, 
unlike in the SM, the {7(1) symmetry of the theory has a physical interpretation in terms 
of baryon and lepton number. Also, the observed parity-odd nature of the low-energy 
weak interactions is seen to be the result of a spontaneously broken symmetry in the LR 
Model. Furthermore, it has been pointed out that, since parity is an exact symmetry of the 
theory at high energies, the strong CP problem has a natural solution in the context of LR 
models a [io]. In addition, the choice of a Higgs triplet to break left-right symmetry at a 
high energy scale leads naturally to a seesaw relation to explain the smallness of neutrino 


masses 


Cl 0. 


Finally, the extended Higgs sector of the theory gives the possibility, at 
least in principle, of breaking CP spontaneously. 

Many authors have focused on the possibility of spontaneous CP violation in LR models, 
both in minimal versions as well as in extended versions containing more complicated Higgs 
sectors 0000 00. The results in minimal versions of the LR Model have 
by and large been somewhat negative regarding the possibility of generating spontaneous 
CP violation, while investigations of more complicated models have had more encouraging 
results. A further complication in the case of the minimal model concerns the presence 
of neutral Higgs bosons having flavour-violating couplings to quarks (Flavour Changing 
Neutral Higgs (FCNH) bosons). Recently, Barenboim et al studied the decoupling limit 
of the minimal CP-invariant Higgs potential in the LR model Q. They found that non- 
negligible CP violation in the vacuum state is generally accompanied by weak-scale non¬ 
standard Higgs bosons, which are ruled out phenomenologically. For CP-conserving vacuum 
states they found that it was possible to produce an acceptable physical Higgs spectrum, 
but only at the cost of extreme fine-tuning. 

Another well-known complication in the Higgs sector of LR models concerns the vacuum 
expectation value (VEV) of the neutral component of the left-handed Higgs triplet, (S R ) ~ 
VlU 6l . Assuming that all dimensionless coefficients in the Higgs potential are of order unity, 
Vl is related to its right-handed counterpart v R through a seesaw relation, vl k 2 /v R , 
where A: is a dimensionful quantity of order of the weak scale. The complication resides 
in the fact that, in general, the left-handed neutrino mass matrix contains a term that is 
proportional to v R . Barring extreme fine-tunings, this term in the mass matrix forces v R 
to be of order a few eV or less, thereby forcing the right-handed scale, v R , to be extremely 
large. To allow for the possibility of an observable right-handed scale, many authors assume 
that Vl = 0, a scenario that can be arranged by disallowing certain terms from the Higgs 
potential. The offending terms in the Higgs potential have the form Tr^0A^0f and may 
be forbidden, for example, by imposing a symmetry under the discrete operation A R —■► 
A l and A R —> —A R . Unfortunately, invariance under this symmetry also disallows the 
Majorana Yukawa couplings that give rise to the seesaw relation for neutrino masses 0. 
Another approach is simply to assume that the theory is embedded in a grand unified scheme 
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that somehow disallows the troublesome terms from the Higgs potential but still allows the 
desired Majorana Yukawa terms. This of course does not solve this problem within the LR 
model framework. 

In this paper we analyze the Higgs sector of the minimal LR Model in the case that CP 
is not a manifest symmetry of the unbroken theory, a scenario that has not received much 
attention. In the most general case, the introduction of explicit CP violation corresponds to 
the introduction of only one complex coefficient into the Higgs potential, i.e. only one extra 
degree of freedom. As we shall see below, the phase associated with this complex coefficient 
significantly alters the dynamics of the Higgs sector, and we find that the presence of CP 
violation in the vacuum state no longer precludes the existence of an acceptable Higgs 
spectrum. 

In our calculation, in order to suppress Vl to phenomenologically acceptable levels, we 
further assume the presence of an extra horizontal symmetry that is broken by a small 


parameter e# 


le t rie p res 

[HEo-E 


An appropriate choice of charge assignments under this symmetry 


suppresses certain coefficients in the Higgs potential and brings Vl down to the eV scale for 
moderately light right-handed scales (vr ~ 15 TeV, say). The model still does contain 
a certain amount of fine-tuning. In particular, one dimensionless coefficient in the Higgs 
potential needs to be “unnaturally” small (of order e/vl). Since we restrict our attention 
to right-handed scales that are only “moderately” large, this fine-tuning is not as unnatural 
as it might be in general. 

Our paper is organized as follows. In Sec. [HI we describe the model and perform the 
minimization of the Higgs potential to obtain relations among the various Higgs VEVs. Sec¬ 
tion lllll contains several approximate expressions for the Higgs boson masses. Section II V I 
contains our numerical results and in Sec. m we offer some concluding remarks. The Ap¬ 
pendix describes the derivations of the approximate expressions for the Higgs boson masses. 


II. THE MODEL 


A. Minimization of the Higgs Potential 


We consider a LR model with a single bidoublet held 0 ~ (2, 2, 0) and triplet fields 
A l ~ (3,1,2) and A r ~ (1,3,2). Under SU(2)l x SU(2)r the fermion fields in the theory 
transform as 'I > L R —»• Ul,r^l,r and the Higgs fields transform as 0 —> Ul4>Ur and A l,r — > 
Ul,rAl,rU' l r. We also define the held 0 = 720 *t 2 , which transforms in the same way as 0. 
The full Lagrangian of the theory (denoted below as LR Model) is invariant under the LR 
symmetry operation 


0^0 t , 0 0 f j A L ^ A R . 


( 1 ) 


The most general LR symmetric Higgs potential that can be constructed from 0,0, A^ and 
A r may be parameterized as follows jl6{, 


V = 


-lA Tr(0 t 


hi 


Tr 


+Ai 


Tr 


n 2 


+ A 2 


Tr 


^ + Tr(ft 

n 2 


+ 


Tr(0 f 


hi 

n 2 


+ Tr(A H Aj,) 


+ A3 


+A 4 {Tr(00 t ) [Trffi) + Tr(0 f 0)]} + pi [Tr (a l aJ,)] + [Tr(A R A^) 


Tr 

2 


*) Trf0 f 


3 

























+p 2 Tr(A L A L )Tr(A[Ay + Tr(A R A R ) Tr(A^A^) + p 3 Tr(A L A^) TrjA^A 
+p 4 [Tr(A L A L ) Tr(A^A^) + Tr(A^A^) Tr(A R A R ) 

+«i{Tr(#t) [Tr(A L Al) + Tr(A R Ay]} 

+«2 { e i52 [Tr^) Tr(A^A^) + Tr^f) Tr(A L A^)J + h.c.} 

+a 3 [Tr^^AiA^) + Tr(0 f 0A R A^) + fa Tr^0A jR 0 t A^ + Tr(0 T A L 0A^ 

+& [Tr(^A R 0t A y + Tr(0tA L 0Ay] + & [Tr^A^A^) + Tr^A^A^)] . (2) 


All of the coefficients in the above expression are real, with the exception of the term 
proportional to 02 , which has a phase $ 2 - ( The coefficient a 2 itself will be taken to be real 
and non-negative in this work.) Papers on this subject have nearly universally set sin 82 
to zero. As pointed out in Ref. [1£], however, setting sin 62 to zero effectively means the 
vacuum state will not break CP. In contrast to this, we will see that it is quite natural to 
have large CP-violating phases in the Higgs VEVs in the case sin<f 2 7 ^ 0. 

We adopt the following notation for the various components of </>, A L and A R , 


0 = 




Al,r — 


L,R / V^2 ^ ; f{ \ 

$L,R ~<5 L,R /V/2 J 


( 3 ) 


The lepton Yukawa couplings that are consistent with the gauge and parity symmetries 
discussed above are 


- ^Yukawa = ^' iL ^ jR + \ F ij (i’'iL Cr 2 A L ^' jL + ^ r Ct 2 ArI^'j^ + h.C. , (4) 

where F, G and H are 3x3 Yukawa matrices, C = ? 7 2 7 ° is the charge conjugation matrix 
and primes denote gauge eigenstates. Upon spontaneous symmetry breaking, the neutral 
components of the Higgs boson fields obtain VEVs, 


(</>) 


h/V 2 o \ , v ( 0 o \ 

0 k 2 e ia /\/2 J ’ ( ^ v L e idL /V2 0 J ’ ( R ' 


0 0\ 

V R /V 2 0 ) ’ 


( 5 ) 


where k \, k 2 , Vl and v R all refer to the magnitudes of the respective quantities. Phenomeno¬ 
logical considerations lead to the conclusion that v R -C ki,k 2 -C v R . Furthermore, ki and 
fc 2 are at the weak scale, fc 2 + k% — (246 GeV) 2 . Gauge rotations have been used in the 
above expressions to eliminate possible phases in the ki and v R terms 0 The four complex 
neutral fields may then be expanded in terms of eight real fields as follows, 

4>\ = «-H 0 ? + fc!)/V 2 , ( 6 ) 

41 = (7 + i<$ + h)e‘°/s/2, (7) 

7 = (7 + >7 + v L )e‘^/V2, ( 8 ) 

7 = (7 + *7 + 7 / 72 . ( 9 ) 

Minimizing the potential in Eq. ((21) with respect to 05 r , <$£, (f)® 1 , 8 R , and S R leads to 
six equations that need to be satisfied by the coefficients of the potential. By manipulating 
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these expressions we obtain expressions for the /i|, 


lA ai (, v 2 l \ o 3 £ 2 

9 - 1 ' 2 J 2(1-a 


4 


2 

+ 


vi 


1 R 

@2 cos e L - /3 3 £ 2 cos(0 L - 2a)] , 

J l-( 


1 + ArJ + Ai (l + £ 2 ) + 2A 4 £' 


( 10 ) 


A _ «2 

vf { 2 cos a 


+ 


+ 


“**+« +cos( “ - «|J+ 4(i-ac 05a l 1 +1) 

r 1 1 e 2 

2A 2 £ cos(2a) + A 3 £ + -A 4 (l + £ 2 ) cos a - 

2 v ' J cos a 

(3i(l - £ 2 ) cos (0 L - a) - 2/? 2 £ cos 0 L + 2/%£ cos (6 L - 2a)] - 

L J 4(1 — t, z ) cos a 


,(H) 


iAs 


v 2 n 


= P i 


i + -Aj^J + 2 ai ( x + £ 2 ) + 


+2o 2 


n 2 

cos(a + 5 2 ) — cos(a — <5 2 ) — 

V R 


1 - vl/vl ’ 


( 12 ) 


where we have defined £ = k 2 jk\ and e = ki/vR. Both £ and e are small quantities in our 
analysis. We also obtain three other equalities that must be satisfied, 

8a 2 £e 2 sin a sin 5 2 


ocxqc; c sin a sii 
(2pi “ p3) — 


vl_ 

Vr 


/?i£ cos (0 L - a) + (3 2 cos d L + /3 3 £ 2 cos(0 L - 2 a) 


eA 


(13) 


0 = /?i£ sin(6> L - o) + (3 2 sin 0 L + /? 3 £ 2 sin(0 L - 2a) , 


(14) 



£ 2 ) (i-uiM) sin ^ 

= {2£sin(6> L - a)((3 2 + /? 3 ) + sin d L + £ 2 sin(6> L - 2a) Pi}^- 
+£ sin a [a 3 (l + vl/vj^ + (4A 3 - 8A 2 ) (l - £ 2 ) t 


vr 


(15) 


Taking 5 2 —> 0 in Eqs. dT31) - (H51) gives Eqs. (A5)-(A7) in Ref. 0, Equation (THfll yields the 
well-known seesaw relation for vl, vl ~ (some couplings) x e 2 VR. Interestingly, p 2 and p 4 
do not appear in the above expressions and are thus not constrained by the first-derivative 
conditions. (The first of these does, however, appear in an approximate expression for one 
of the Higgs boson masses and is thus constrained through a second-derivative condition.) 

A very important feature of the above equations is that the left-hand side of Eq. m is 
generically non-zero. This gives us considerably more freedom when solving the minimization 
conditions compared to the situation considered throughout much of Ref. jl8| (in which it was 
mostly assumed that sin S 2 = 0). In particular, in our case we can easily obtain a 3 = 0(1) 
instead of 0(e 2 ). This observation is crucial, since the neutral Higgs bosons with flavour non- 

diagonal couplings to quarks have masses ~ yja 3 v 2 R /2 (see Eq. (|2Hj) below). If o 3 = 0(1), 
these FCNH bosons attain masses at the scale vr and are phenomenologically viable. If, 
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however, a 3 is of order e 2 , then the FCNH bosons have masses at the weak scale, a scenario 
that presents considerable phenomenological problems. Let us examine this argument a bit 
more quantitatively. Suppose for the moment that sin S 2 = 0 and Vr/vr -C 1. Suppose 
furthermore that sin a 7 ^ 0, so that the vacuum state breaks CP. In order to satisfy Eq. (11511 
in this case, one must have a 3 ~ ( 8 A 2 — 4A 3 )(1 — £ 2 )e 2 ; i.e., unless A 2 is very large, a 3 must 
be of order e 2 , leading to phenomenological difficulties. If sinh 2 7 ^ 0, however, things are 
quite different. I 11 this case Eq. m becomes 


2 a 2 sin 5 2 ~ a 3 £ sin a 


(16) 


(where we have assumed that Vr/vr is negligibly small and that the various Higgs coefficients 
are not anomalously large). It is possible to satisfy this expression with a 3 = 0(1). In 
particular, we can use the above expression to determine the phase of (4> 2 ), 


a ~ sin 


-1 


, 2 ck 2 sin<5 2 ' 

. « 3 ? . 


(17) 


With the choice we make below for the charges under the approximate horizontal symmetry, 
we have a 3 ~ 0(1) and a 2 ~ 0(£). In this case a is of order unity if S 2 is. To summarize, if 
sin S 2 is non-zero, one can in fact have a vacuum state that violates CP while simultaneously 
evading the FCNH problem. This possibility was noted in Ref. |18], but was not explored 
in detail. 


B. A Broken Horizontal Symmetry 

Let us return to the seesaw relation for the VEVs derived from Eq. m and noted above. 
If all the coefficients in the Higgs potential are of order unity one has Vr ~ (some couplings) x 
6 2 vr (note that 6 2 vr = kf/vR). If k\ is of order the weak scale and vr is of order 10’s of 
TeV, then vr is of order a few GeV. Such a value for vr is extraordinarily large from the 
point of view of neutrino physics, an observation that can be understood by examining the 
3x3 mass matrix for the light, mostly left-handed neutrinos. This mass matrix is given by 
the following approximate expression, 

M y ~ M RL — MrrM rr MJ) Ri ( 18 ) 

where 

Mrr = (Gki + Hk 2 e~ ia )/V2 , M LL = Fv L e i(>L /V2 , M RR = Fv R /V 2, (19) 

with F, G and FI being 3x3 Yukawa matrices (see Eq. O). Phenomenologically, the 
elements of M u must be at the eV scale (or smaller). The second term in Eq. (fTHl) is the 
usual “seesaw” term - it scales roughly as kf/vR, thereby suppressing the neutrinos’ masses. 
If wr is only “moderately” large (of order 10’s of TeV), the Yukawa matrices G and H need 
to be suppressed in order that the elements in this term are at the eV scale. 1 The first term 
in Eq. m also benefits from a seesaw relation in that it is proportional to vr ~ kl/v R . The 


1 Such a suppression can be obtained in the horizontal symmetry scheme that we adopt in this paper |2 
although we do not consider neutrino masses further here. 
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seesaw suppression is not enough if vr = 0(10 TeV), however, since then vr is of order a 
few GeV (as noted above). To solve the problem, one could suppress Mll by forcing the 
elements of F to be very tiny (of order 10~ 10 , say), but then the right-handed Majorana mass 
terms ( Mrr ) would also become very small (of order keV for vr of order 10’s of TeV). This 
would render the seesaw mechanism for neutrino masses inoperative (in fact, Eq. m would 
become invalid) and would generically yield GeV-scale neutrinos. To remedy the situation 
one needs to force the coefficients fa, in Eq. m to be very small or zero, thereby forcing vr 
to be very small or zero. The fa can be made to disappear via an exact symmetry, but, as 
shown in Ref. 0 . the disappearance of the fa comes at the expense of Majorana or Dirac 
mass terms in this case. 

A simple approach that suppresses the fa without eliminating them completely is to adopt 
an approximate horizontal 17(1) symmetry that is broken by a small parameter e# 


o adopt 

BE, 


21]. In this approach, each of the fields in the theory is assigned a charge under the symmetry, 


and terms in the Lagrangian that are not singlets under this symmetry become suppressed 
by €fj, where n is some combination of the charges of the fields involved. We adopt the 
charge assignments 


21 


Q(Ai) = -Q(A„) = 

Q(4>) = - 2 , 


( 20 ) 

( 21 ) 


with the result that the dimensionless coefficients in the Higgs potential scale as follows 


ai,a 3 ~0(l), a 2 
Pi~0{e%), fa 
A 1} A 3 ~ 0(1), A 2-0(6%) 

P1,P2,P3 ~ 0(1), 


o(0 H ) ~ 0(0 

0($), fa~0(e 

0(0), A 4 ~ 0(e 

/?4 ~ 0(e 




2 ) 


0(0 


( 22 ) 


In the numerical work below we set £r = 0.3 and £ = k 2 /k\ = 3/181 ~ 0.017 0 , so 
that ejj ~ £, as noted in the expressions above. Furthermore, for a moderate right-handed 
scale (of order 15 TeV), e = O/vr ~ (246 GeV/(1.5 x 10 4 GeV) ~ 0.016, so that e ~ £. 
(Of course, for very large right-handed scales e becomes much less than £.) The fa are 
particularly suppressed in this scheme, with the result that the VEV seesaw relation now 
becomes 


v L ~ e™e\ R . 


(23) 


Taking vr = 0(15 TeV), we have vl ~ 0.1 eV, which is phenomenologically acceptable. 


C. Approximate Expressions for the Minimization Equations 

It is useful to consider approximate versions of the six minimization equations (Eqs. (II011 - 
m which can be obtained by expanding in the small parameters £, e and vr/vr (as well 
as taking into account the charge assignments for the various Higgs potential coefficients). 
Of these, we consider e and £ to be of approximately the same order (approximately 0.017), 
while vr/vr is miniscule (of order 10 -14 ). We also bear in mind the scaling of the various 
coefficients in (12211 . recalling that e 4 H ~ £. An approximate (and slightly rearranged) version 
of Eq. (fT51) has already been given above in Eq. (116(1 . (The corrections to this expression 
are of order e 2 £ and £ 3 .) Furthermore, Eq. (ED does not need to be expanded at all, since 
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each of the three terms is of approximately the same order. From the other four equations 
we obtain the following expressions, 


ai 

v 2 r 

i4 

v 2 r 

hi 


J R 


y-i« 3 e 2 + A y + oty.ev), 

a 2 cos (a + <5 2 ) + + £e 2 ) 


cos a 


4 cos a 


v L = 


Pl +-a l£ 2 + (9(£ 2 e 2 ), 

[/?i£ cos(6U - a) + (3 2 cos 9 l + /C^ 2 cos(6U - 2a)] e 2 u/j 

(2pi - Ps) 


1 + 0(£V) 


(24) 

(25) 

(26) 
(27) 


Let us hrst consider the decoupling limit for the Higgs fields; i.e. let e become very small, 
keeping £ and e# fixed. Suppose that we are given a particular model for the Higgs potential, 
so that all the coefficients in the Higgs potential are fixed. In the decoupling limit, Eqs. m 
and (1231-0 represent four equations in three unknowns (vr, £ and a). The system is 
overconstrained and there must be fine tuning (at order e 2 ) among the coefficients in the 
Higgs potential. This fine-tuning is quite severe if Vr is (R(10 7 GeV) (as is often assumed in 
the decoupling limit), in which case e 2 = 0(1O -9 ). The remaining two equations (Eqs. (fTHl 
and (|27|) 1 provide constraints on the remaining three parameters (ul, e and 6l)- 

It is well-known that the Higgs sector of the LR Model contains fine-tuning. In some cases, 
this fine-tuning extends to more than one equation that must be satisfied among coefficients 
in the Higgs potential 0- Having said this, we note that the fine-tuning problem becomes 
less severe if vr is taken to be at a moderate scale; i.e., of order 15 TeV. Suppose again 
that one is given a set of coefficients in the Higgs potential and that one wishes to solve the 
approximate expressions in Eqs. (ED, ED and ED-ED for the six unknowns Vr, Vl, Or, e, 
£ and a. One could proceed as follows 2 : 


1. Use Eq. (I2H1) (with e ~ 0) to determine vr. 

2. Use Eqs. ED and ED to determine a and £ = k 2 jk\. 

3. Use Eq. ED to determine 9r. 

4. Use Eq. ED to determine e = ki/vR (thereby determining the weak scale). 

5. Use Eq. ED to determine Vr. 

The fine-tuning issue surfaces when one uses Eq. ED to determine e = ki/vR. If the 
dimensionful parameter p 2 is taken to be “naturally” of 0(v R ) (as is pg, from Eq. ED), 
then ai must equal 2 P \/v R to 0(e 2 , £ 2 ) ~ 10“ 4 ; that is, one requires an 0(e 2 , £ 2 ) cancellation 
between two quantities that are naturally of order unity. Alternatively, one could suppress 
p 2 “by hand,” so that it is of order k 2 = e 2 v 2 R . In this case ct\ also needs to be suppressed 
by hand so that it is of order e 2 . In our numerical work we adopt the latter approach 3 . 


2 This procedure could be made more precise by using the full expressions for the six first-derivative con¬ 
straint equations. 

3 At first glance it might appear that another source of fine-tuning occurs in Eq. (l25l) . since the terms on 

the right-hand side of the expression are of order £. Such is not the case, however, since scales like 

e R ~ £, so a “natural” choice for is /i| ~ ^v R . 








III. HIGGS BOSON MASSES 


In this section we give approximate expressions for the masses of the various Higgs bosons. 
These expressions help determine phenomenologically viable scenarios for the Higgs sector 
of the LR Model. The Higgs boson mass matrices are detemined by taking double partial 
derivatives of the Higgs potential with respect to the relevant Higgs fields. The neutral mass 
matrix is 8 x 8, real and symmetric. The singly- and doubly-charged mass matrices are 
both Hermitian and are 4x4 and 2x2, respectively. The eigenvalues of these matrices 
correspond to the squares of the various Higgs boson masses and they must be positive. Both 
the neutral and the singly-charged mass matrices have two zero eigenvalues, corresponding 
to the Goldstone modes absorbed to give masses to the Z and W bosons, respectively. 

Equations below give approximate expressions for the squares of the Higgs boson 

masses. The derivation of these expressions is explained in more detail in the Appendix. 
For the neutral Higgs bosons we obtain the following approximate expressions (aside from 
the two massless Goldstone modes), 


m 

m 


(0)2 

SM 

(0)2 
OL 3 


< n 

mjj ’ 2 


2A 1 (1-K 2 )+8A 3 £ 2 

+ 8A 4 £ cos a kf , 

(28) 

«3(l + £ 2 ) 2 

[2 (i-a J 

v 2 r (two) , 

(29) 

~ Pi + 4a 2 e 2 £ sin a sin h 2 j Vr (two) , 

(30) 

2pi4 , 


(31) 


where the word “two” in parentheses indicates pairs of degenerate or nearly degenerate 
eigenvalues. Note that two of the neutral mass-squared eigenvalues do depend somewhat on 
a, the CP-odd phase that appears in (0) (see Eq. (I5jl ). Explicit reference to the phase h 2 
(the sole CP-odd phase in the Higgs potential) may be eliminated in the expression for m® 2 
by using the first-derivative condition in Eq. m We have left the expression as is since it 
is more compact. For the singly-charged Higgs bosons we obtain the following expressions 
for the two non-Goldstone modes, 


rriM 2 ~ 
U a 3 ~ 


«3 

T 


(t-K 2 ) ^ 

( 1-0 2 1 * 


J R > 


m 


IP -m<r + ^ 2 (l-yw, 


and for the doubly-charged Higgs bosons we have, 


m 


<r 12 - m<f +f e 2 (i - ey R , 


m 


P3 

(++)2 

P2 


2 P2 + ye 2 (l - e) 


VR ■ 


(32) 

(33) 


(34) 

(35) 


The expected corrections to the above expressions are described in detail in the Appendix. 
At this point we simply note that for the range of parameters considered in this work 
the expressions for the singly- and doubly-charged Higgs bosons are essentially exact, as is 
Eq. m l° r two °f ll ie neutral Higgs bosons. The other expressions for the neutral Higgs 
bosons have varying levels of accuracy. Reference [17] contains approximate expressions for 


9 












the Higgs boson masses assuming vr = 0 and a — <5 2 — 0. The above expressions agree with 
those in Ref. 0 in the stated limit, and neglecting some higher-order corrections. Also, to 
a very good approximation, rri^J 2 + m^ + - )2 = 2m^ 2 . 4 

A few comments are in order. First of all, the neutral Higgs boson with essentially SM- 
like, flavour-diagonal couplings to the quarks (except for small corrections) also generically 
has a mass that is at the weak scale. Furthermore, as was noted above, the six first-derivative 
equations allow a 3 to be of order unity, so the masses of the FCNH bosons (Eq. (I2H1) - ) are 
comfortably of order Vr. In fact, in the horizontal symmetry scheme that we employ, Ai, 
a 3 , Pi, P 2 and p 3 are all of order unity, so all non-SM Higgs bosons have masses that are 
generically of order vr. This situation is to be contrasted with the case in which the Higgs 
potential is CP-invariant (sin <f 2 = 0). I 11 that case it is difficult to maintain a separation in 
scales between the SM-like Higgs boson and the non-SM-like Higgs bosons [18]. 

In addition to enforcing the first-derivative conditions (Eqs. (fTUlh (FT51P . we also need to 
ensure that the extrema of the potential are actually minima. Physically, this reduces to the 
requirement that all of the mass-squared eigenvalues be positive. To a good approximation 
it is sufficient to require that Ai, a 3 , p± and p 2 all be positive and that p 3 be larger than 2 p 1 . 
The value of the potential at the minimum is approximately — piuJj/4, which is negative, 
since p\ is positive. 


IV. NUMERICAL RESULTS 


In this section we perform a numerical analysis of the model to determine if there are in 
fact combinations of coefficients in the Higgs potential that give phenomenologically viable 
spectra of Higgs bosons. We also wish to determine whether there might be correlations 
among the magnitudes and phases of the various VEVs, or whether there might be con¬ 
straints on the sizes of the phases a and 0 L . Such correlations or constraints would affect 
quark and lepton mixings and could have important phenomenological consequences in K-K 
and B-B mixing and in neutrino physics. 

Our basic approach is to select sets of coefficients consistent with Eqs. m-m and in 
the ranges specified by Eq. (j22l) (except for aq, which is taken to be of order e 2 , as noted 
in Sec. m and then to compute the masses of the various Higgs bosons by diagonalizing 
the mass matrices numerically. This procedure also allows us to check the reliability of the 
approximate expressions for the masses in Eqs. ©-©• 

As noted above, the mass-squared eigenvalues for physical Higgs bosons must be positive, 
which essentially forces Ai, a 3 , pi and p 2 all to be positive and p 3 to be larger than 2 p 1 . 
Let us consider some further constraints that follow from direct and indirect limits on the 
masses of the Higgs bosons. 

The SM-like Higgs boson has couplings to quarks that are very similar to those of the 
Higgs boson in the SM, so we may use limits on the SM Higgs boson to constrain Ai. The 
Particle Data Group suggests the range 114.4 GeV< m® r < 193 GeV for the SM Higgs 
boson based on a global fit 0, which translates into the range 0.1 1$ Ai ;$ 0.3. Direct lower 


A similar relation was reported in Ref. [2l| . although our expressions differ in some respects from the 
ones contained in that paper. For example, the authors of Ref. 1 2 .‘ll | include an extra term in the Higgs 
potential, and their expressions for the neutral mass matrices appear to be incorrect by an overall factor 
of 2. See also Ref. Q- 
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TABLE I: Ranges used in the numerical analysis for the dimensionless coefficients in the Higgs 
potential, assuming vr = 15 TeV. The horizontal symmetry breaking parameter eu is set to 0.3 in 
our numerical work. All coefficients not included in this table scale as in Eq. & with minimum 
and maximum values of =f 2 multiplied by eu to the appropriate power. Thus, for example, we 
require that — 2e^ < (3\ < 2e^. 


coefficient 

minimum value 

maximum value 

V 

0.1 

0.3 

a 3 

0.2 

2 

Pli P2 

0.06 

2 

P3 

2pi + 0.2 

2 

Oil/ e 2 

-2 

2 

a 2 /e A H 

0 

2 

52 

0 

2tt 


bounds from LEP experiments on the singly-charged and doubly-charged Higgs bosons are 
typically of order 100 GeV 01 0)0)0]; while indirect searches can produce a much higher 
reach. For example, under certain assumptions regarding Yukawa coupling constants the 


igconsta 

lower bound on doubly-charged Higgs boson masses extends to the TeV range [0,0]. Future 
indirect probes of the TeV-range masses are warranted in other high and low-energy (Mpller- 
scattering) experiments. In the neutral case, the Higgs bosons with masses determined 
primarily by a 3 have flavour non-diagonal couplings to quarks and can contribute to neutral 
meson mixing at tree level. Limits from such processes are somewhat model dependent, but 
can provide lower bounds in the range of several TeV and even up into the range of tens 
of TeV 0, HH 0 , 0]. For the purpose of our numerical analysis we fix vr = 15 TeV (as 
well as k\ ~ 246 GeV and £ = k 2 /k\ — 0.0166) and place lower bounds on a^, pi, P2 and 
p 3 in such a way that all non-SM Higgs bosons obtain masses of order 5 TeV or higher. For 
the upper bounds on these coefficients we somewhat arbitrarily choose the value “2,” which 
is consistent with “order unity.” The bounds employed for all coefficients are described 
in Tabled The only coefficient that does not scale as in Eq. (1221 ) is oq, which we have 
suppressed “by hand” as discussed in Sec. Ill Cl 

Our specific procedure for obtaining VEVs and sets of coefficients that satisfy the six 
constraints in Eqs. (CDD-dlSD is as follows. First, we choose a random value for vr (between 
zero and 2 ef/kf/vR), as well as random values for the phases a and Or (in the range zero 
to 2n). In a similar manner we choose random values for S 2 and all other coefficients in 
the Higgs potential except for a 2 , p 3 , the pf and one of the The six first-derivative 
equations (Eqs. (fTUl) - (fT51) ) are then used to compute these remaining six coefficients. If a 2 , 
p 3 and the one remaining /3{ are all in acceptable ranges (as determined by Table |T|) we 
proceed to compute the various mass-squared eigenvalues. We do not enforce any particular 
constraints for the p 2 . Numerically, p 2 is positive and of order v 2 Rl as one would expect from 
Eq. (I2H1) (recall that p\ is positive). p\ is numerically of order the weak scale squared in 
magnitude and p\ is intermediate between the squares of the weak and right-handed scales 
in magnitude. p\ and p% can be positive or negative, as one might expect from Eqs. (1241) and 
(1251) . This result may seem counterintuitive compared to the SM. In the LR model under 
consideration it is unusual for the Higgs potential to be a local maximum at the origin. 
The origin is generically a generalized saddle point in the sense that the eigenvalues of the 
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FIG. 1: Frequency plot for the masses of the lightest Higgs bosons. The upper left plot shows 
the results for the SM-like Higgs boson, while the other plots show results for the lightest neutral, 
singly-charged and doubly-charged nonstandard Higgs bosons. 


second-derivative matrix generically have mixed signs. In any case, the Higgs potential is 
never a local minimum at the origin, since /i| > 0. 5 

Figures |U and [21 show the results of our numerical analysis. The first of these shows a 
frequency plot of the lightest Higgs boson masses. While we have not imposed any direct 
mass constraints on the numerical results, the ranges described in Table |T] effectively limit 
the mass of the SM-like Higgs boson to be in the range 100 GeV < mgjq < 200 GeV and 
those of all nonstandard Higgs bosons to be of order 4.7 TeV or greater. The scatter plots 
in Fig.[2]show various correlations between vl , 0 l, ol and 82- The correlation between a and 
62 is easily understood by examining the approximate relation in Eq. (USD - since « 2 , ct 3 and 
£ are all positive quantities, sin 82 and sin a must have the same sign. The remaining three 
plots show various combinations of a, vl and 6 r, quantities that affect quark and lepton 
masses and mixings. In particular, notice that the phases a = arg ((</>!])) an d @l = arg ((5°)) 
are not constrained to be small. Also, vr is of order 0.1 eV, which is phenomenologically 
viable. From these plots it does not appear that there are strong correlations between a, vr 
and 6 r. A numerical investigation of the case vr = 50 TeV yields similar plots (except that 
vl is reduced in magnitude because it scales as 1 /vr). 

Our numerical analysis also allows us to check the reliability of the approximate expres- 


5 In the eight-dimensional space spanned by the neutral fields bib etc. (with all VEVs and phases set to 
zero), a local maximum at the origin is obtained if > 0, 2y| < /jl\ and — 2/z| < Mi- This situation is 

“unusual,” since is typically much smaller in magnitude than is /A ■ The condition for a local minimum 
at the origin is obtained by reversing the inequalities in the three expressions. 
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FIG. 2: Various correlations between a, 62 , 0l and vl- 


sions for the mass-squared eigenvalues in Eqs. (EBP-(ESI)- The eigenvalues were computed 
numerically from the mass matrices and compared to the approximate expressions. For the 
singly- and doubly-charged Higgs bosons, the “exact” (numerical) values agreed with the 
approximate expressions to within a few parts in 10 15 for vr = 15 TeV. The level of agree¬ 
ment for m® 2 was similar. Corrections for the other four neutral mass-squared eigenvalue 
expressions are described in the Appendix. 


V. DISCUSSION AND CONCLUSIONS 

We have considered the Higgs sector of the minimal Left-Right model with explicit CP- 
violation in the Higgs potential in the “moderate” decoupling limit, i.e., when the scale set 
by the VEV of the right-handed Higgs field vr is in the range 15-50 TeV. This intermediate 
regime provides (at least in principle) testable effects of the remnants of RH symmetries in 
upcoming collider experiments. At the same time, inclusion of explicit CP violation in the 
Higgs potential allows for the generation of a viable spectrum of Higgs bosons containing 
one SM-likc Higgs boson with mass of order the weak scale and several heavy neutral and 
charged Higgs bosons with masses of order vr. Supplemented by an additional [/(1) horizon¬ 
tal symmetry, this model alleviates most of the fine-tuning effects associated with minimal 
LR models in the decoupling regime 0 , resulting in an SM-like low energy effective theory. 
Yet, some small amount of fine-tuning is still required. We have also performed numerical 
simulations of the Higgs spectrum in our model, confirming our phenomenological analysis 
and the reliability of the adopted approximations. The power counting for the Higgs poten¬ 
tial coupling constants provided by the broken 17(1) charge assignments also allows for the 
natural generation of neutrino masses in such a model. 
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APPENDIX: GOLDSTONE BOSONS AND APPROXIMATE MASSES 

In this appendix we give exact expressions for the charged and neutral Goldstone bosons. 
These expressions are used to construct matrices that perform exact block diagonalizations 
for the singly-charged and neutral mass matrices, separating out the “zero” eigenvalues 
for the Goldstone modes. We also explain our procedure for obtaining the approximate 
mass-squared eigenvalues given in the text (Eqs. (I2%l) - (135jl ). 


1. Goldstone Bosons 


The kinetic terms for the Higgs boson fields in the Lagrangian are given by 7} 


£id„ = Tr 


+ Tr (D^A^D^Al + Tr (D M A R ) f D»A R 


(A.l) 


where covariant derivatives are defined as 


= B,A - f WW + . 


D „A 




— d^A l,r — 


W_ 

2 L 


w l,r^ t<1 ^l,r -ig'Bfj A l ,r 


(A.2) 

(A.3) 


Inserting Eqs. © and ©-© into Eq. KB yields the mass matrices for the gauge bosons, 
which may be diagonalized to obtain the physical mass eigenstates (Wf^, Z\ )2 and A^) in 
terms of the original gauge bosons associated with SU{2)l x SU(2)r x U(1)b-l (wH = 
(li'X,' T / A, W RR and B^). Equation ilA.lli also yields bilinear couplings between 

the physical gauge bosons and various linear combinations of the Higgs fields. One such term 
is proportional to igWAd^Gf , for example, and allows one to identify Gf as a (would-be) 
Goldstone boson. Proceeding in this manner one obtains expressions for the two charged, 
orthogonalized Goldstone bosons, 


Gt 

Gt 


1 



(ki sin ( — k 2 cos () e *“0^ + (k i cos ( — k 2 sin () 0^ 

- V2v R e- ia sm( <5+ - V2v L e~ i6L cos (S+] , 

(fci cos ( + k 2 sin () e~ ia 4>i + (— k\ sin ( — k 2 cos C) <t> 2 

— \f2v R e~ ia cos (5ft + V2vLe~ ldL sin (5 J , 


(A.4) 


(A.5) 
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where 7V0 2 are normalization constants and 


sin 20 


2 A; i A; 2 

'(l>K - I’D 2 + 4 1/2 


(A.6) 


We may also construct two other normalized, orthogonal linear combinations of the four 
holds 00, 00, 50 and 50, 


<t>t 

<t>t 


1 


1 

K 


(V iQ 00 + 

(r£e ~ ia (j)i 


r$ 2 44 + r+e-^Sft) , 

(A.7) 

+ r0 2 0+ + r+e-* Q 5+ + e-^50) , 

(A.8) 


where the are simple functions of k \, k 2 , v R and Vl- (Forcing 0+ to be orthogonal to G0 
and G0 determines r 32 and r 3?J ; similarly, forcing 00 to be orthogonal to G0, G0 and 00 
determines r0 L , r0 2 and rf 3 .) Putting the above expressions together, we form the unitary 
matrix 77.T) that relates the two bases, 


(Crt\ 


(<Pt\ 

Crt 

00 

= TZ {+) 

<t>2 

& + R 

) 


U L J 


(A.9) 


The matrix H( + ) may be used to bring the 4x4 singly-charged mass matrix into block- 
diagonal form, with a non-zero 2x2 block in the lower-right, 


1Z {+) * M 2 + 1Z {+)T 


0 0 \ 
0 A40 2X2 / 


(A.10) 


One can similarly block-diagonalize the neutral mass matrix, although the diagonalization 
of the neutral gauge bosons is somewhat more involved than that of the charged gauge 
bosons. Furthermore, both g and g' are involved, so that the physical mass eigenstates 
depend on the Weinberg angle (tan 2 9w = g ' 2 / ( g 2 + g' 2 )). Since our purpose is to block- 
diagonalize the neutral Higgs sector, and since the Higgs sector has no dependence on 
9\y, we will consider the (unphysical) limit g' / g —> 0. This means our expressions for 
the Goldstone modes, while useful for diagonalization purposes, will not correspond to the 
actual combinations of Higgs bosons “eaten” by Z\ 2 . Taking this limit, and proceeding as 
in the charged Goldstone case, we obtain 


G? 

G° 


where 


1 


N? 

1 




(cos 0 + sin 0) (ki4>i l ~~ &202*) + 2 v R sin 050 — 2 v R cos 050 , 
(cos 0 — sin 0) (fci00 — /c200) + 2vr cos 050 + 2vl sin 050 , 


(A.ll) 
(A-12) 


sin 20 


kf + k 2 

[4 {v 2 R ~vl) 2 + (kl + kl) 2 } 1/2 


(A.13) 
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Defining orthogonal combinations 


€ = ^ (W? + O , (A-14) 

« = jjjj (r>? + + r»,i“ + Sf) (A.15) 


(where the rZ are again determined by forcing the various fields to be orthogonal), we 
construct the orthogonal matrix TZ^ that relates the two bases, 


( G °l) 


(<ff\ 

g 2 

€ 

= TZ {0) 

<t>2 

W&V 


U?/ 


Defining an 8 x 8 orthogonal matrix 


n {0) 


0 A 

14X4 0 )’ 


we then have 


n w M 2 0 n^ T 


0 0 \ 

0 -^0,6x6 ) 


(A.16) 


(A. 17) 


(A.18) 


where Alp 6x6 is the symmetric 6x6 mass matrix in the basis $o = (0a\ 4>b\ ■> ^ ; 8r, ^l)- 


2. Approximate Mass-Squared Eigenvalues for Physical Higgs Bosons 

Throughout this paper we assume that vr is small - of order 0.1 eV or less - so that 
it is at a natural scale for neutrino masses. In this paper we also assume that the $ are 
suppressed through a horizontal symmetry, as indicated in Eq. m, although any other 
scheme that introduces numerical scaling similar to Eq. m is indeed acceptable. Such a 
suppression of the fp leads naturally to the assumed small value for vr, even if vr is only 
“moderately large” (vr ~ 15 TeV), as indicated in Eqs. (1231) and (1271) . It is important in 
our analysis that vr and the are not identically zero, since vr and its associated phase, 
Or, are significant parameters for neutrino physics (see Eqs. m and (1X21)1. Nevertheless, 
for calculating Higgs boson masses, it is an excellent approximation to consider the limit 
Pi, Vr —> 0. This statement will be quantified in the following. 


a. Doubly-Charged Higgs Bosons 

The doubly-charged Higgs boson masses are determined by a 2 x 2 Hcrmitian matrix. 
The off-diagonal elements in this matrix are proportional to p^VrVr, Pik\k 2 , P 2 k^ and P^k\. 
Assuming the scaling given in Eq. (1221) . the largest of these terms scales approximately 
as 0(e 2 e 1 R v R ). The diagonal elements are of order v 2 R . Assuming no extreme accidental 
degeneracies in the diagonal elements, the corrections to the mass-squared eigenvalues due 
to the off-diagonal contributions are of order e 4 e 2 R v R . With e ~ 0.016 (for vr ~ 15 TeV) and 
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r••u 


6h = 0.3, ~ 0(10 20 ), and one can safely take the limit /3j,V£ —> 0 for the off-diagonal 

elements. Taking this limit for the diagonal elements is also a very good approximation, since 


it only introduces a tiny error (of order e 


4 _40„,2 




V R ~ 


10 28 X-Vft). The approximate eigenvalues 


are thus simply the diagonal elements of the doubly-charged mass matrix in the limit fii,VL —> 
0. These approximate eigenvalues are denoted by m^ + ^ 2 and in Eqs. (USD and USD 

and correspond to the fields 5~£ + and 5^ + , respectively. Our numerical work yields excellent 
agreement between these approximate expressions and the values obtained numerically by 
diagonalizing the full mass matrix inl5 

numerical study with vr = 15 TeV. 


the values agree to within a few parts in 10 15 in our 


b. Singly-Charged Higgs Bosons 


For the singly-charged Higgs bosons it is also an excellent approximation to consider the 
limit fliiVi —> 0. The justification for this approximation is somewhat more intricate to 
argue than in the doubly-charged case, but numerically it does appear to be an excellent 
approximation. (The numerical agreement for vr = 15 TeV is similar to that found in the 
doubly-charged case.) Taking the limit f3{, vl —> 0 in the mass matrix and in TZ >k+ \ one finds 
that 5r decouples from the other three singly-charged Higgs fields, yielding 


/ 


Mi ~ 


m 


m 2 + S,e ia 

T2« l - ^ 
0 


ia -^m 2 + e(l — £ 2 ) 0 N 

m+£ 2 ^m+e£(l-£ 2 )e io 0 

>^((l-( 2 )e-” |m 2 e 2 (l-e 2 ) 2 0 

0 0 m (+) 2 ) 


(A.19) 


in the basis ((ftf, $r, $l ), where 


m 


2 

+ 


Vr 

2 (i-a 


(A.20) 


and where the approximate eigenvalue trM 2 is given in Eq. m- Performing the unitary 
transformation in Eq. (IA. 101) yields the approximate expression for the remaining non-zero 
eigenvalue, denoted m^ 2 in Eq. This latter eigenvalue is associated with the held 


_i_ ca+ + e (i _ ^)e~ ia 5^/V2 /N+. 


c. Neutral Higgs Bosons 

Calculation of the approximate neutral mass-squared eigenvalues is simplified by once 
again taking the limit (3i,VL —■► 0. We begin with the orthogonal rotation shown in Eq. (| A. 1 81) 
(taking the limit Vl —> 0 in 1ZM). This rotation yields the 6x6 matrix Mq 6xS , whose 
eigenvalues correspond to the six non-Goldstone Higgs bosons. In the limit /?*, vl —■> 0, two of 
the fields (<5° r and ~ 5°^) decouple from the rest, allowing for immediate identification of 
their mass-squared eigenvalues. The eigenvalues are degenerate in this limit and are denoted 
by m^ 2 in Eq. In numerical tests with vr = 15 TeV, the approximate expression for 
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these eigenvalues has an accuracy comparable to those for the singly- and doubly-charged 
masses. Removing the entries corresponding to theseytwo fields reduces the size of the 
remaining mass matrix to 4 x 4. We call this matrix Adg 4x4 in the (0°*, 0 4 r , (ffi, 5 R ) basis, 
where is defined in Eq. (1A.14I) . The procedure so far can be illustrated schematically as 
follows, 


Ml 




-^0,6x6 


remove $r 




A/I 2 

- /Vl 0,4x4 ' 


(A.21) 


The helds and S R are approximately mass eigenstates, with the remaining two approxi¬ 
mate mass eigenstates being the following, 


0SM — 


/Or 


v / r+e 

i 

7t+f 


(K + w°2 r ) , 


Dehning 


= 


/ o 

i 

o 


v^-K 2 

vA+4 2 

0 

0 

f 

i 


o\ 

o 

0 


\0 0 

the mass matrix in the basis (0sm, 4>a\ i $r) i s 

/ 


u 2 — 7?1°) a A 2 
• /vt 0,4x4 — /v 2 - /Vl 0,4x4 


1Z^ T 



0 lj 


is 



e 2 

e 2 Z e 2 £ 

eoq,e£ 2 \ 

e 2 £ 

1 e 2 £ 2 



e 2 £ 2 1 



\eai,e0 e£ e£ 


'R ) 


(A.22) 
(A.23) 


(A.24) 


(A.25) 


where the expression on the right gives the orders of magnitude of each of the entries, taking 
into account the scaling in Eq. (1221b In the 1-4 and 4-1 entries we have explicitly included 
a term of order eaq, which could be of order e if one chooses to take aq ~ 0(1), but would 
be similar to other terms in that entry if one would choose aq rs./ 0(£ 2 ) or 0(e 2 ). In the 
following we assume that aq rv./ 0(e 2 ). 6 

The diagonal elements of -Mq 4x4 provide good estimates of the eigenvalues of the matrix. 
One could in principle compute corrections to these estimates by assuming that Ado 4 x 4 
is diagonalized by a “small” unitary rotation (i.e., by a unitary matrix that is essentially 
unity along the diagonal and that has small off-diagonal elements). Such an approach is 
complicated by the fact that the 2-2 and 3-3 elements in AdQ 4x4 are nearly degenerate. In 
fact, the difference between the two elements is of order e 2 ev 2 R , which is the same order of 
magnitude as the 2-3 element of Adj( 4x4 . This situation can lead to a large amount of mixing 


6 Upon diagonalization, the “an” term in the 1-4 element in the mass matrix makes a contribution of order 
af e 2 Vfi to the SM-like mass-squared eigenvalue m^ 2 . If on ~ 0(e 2 ), this term is negligibly small and may 
be ignored. If a± ~ 0( 1), however, the contribution is similar in magnitude to the leading contribution 
(~ 2Aie 2 ufj. - see Eq. (128fl l and it must be taken into account carefully. 
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in the 2-3 block upon diagonalization of and can also complicate matters somewhat 

for computing corrections to the 1-1 and 4-4 elements. We choose instead simply to use the 
diagonal elements to estimate the eigenvalues. 

The approximate eigenvalue denoted m g^ 2 in Eq. (I2%1) corresponds to the field 0sm in 
Eq. (IA.22D . Corrections to Eq. (IA.22I) are expected to be of order e e v Rl e A ^ 2 v R and e 2 ^ A v 2 Rl 
where we have assumed that the various Higgs potential coefficients scale as in Eq. ( 1221 ) 
(except for on, which is taken to be of order e 2 ). The approximate eigenvalue denoted m ® 2 
in Eq. dD corresponds to the nearly-degenerate fields and (j )® r . Corrections to Eq. ©) 
are expected to be of order Finally, the approximate eigenvalue corresponding to 

5 r is denoted by m ® 2 in Eq. (EH). Corrections to this expression are also expected to be 
of order e 2 ev 2 R . In our numerical study with vr = 15 TeV we have compared the exact and 
approximate expressions for the various mass-squared eigenvalues. The numerical differences 
are consistent with the quoted corrections for mgy 2 , m ® 2 and m® 2 , although sometimes 
the corrections are larger due to accidental degeneracies or combinations of coefficients that 
occassionally give large enhancements. An example of the former effect occurs if the 2-2 and 
4-4 or 3-3 and 4-4 elements of M.q 4x4 are nearly degenerate. In that case, corrections to the 
respective eigenvalues can be of order e^v 2 R instead of e 2 ev 2 R . An example of the latter effect 
occurs for m^ 2 , for which one of the leading corrections is proportional to 32 X 2 ^ 2 4/«3- 
If A 3 ~ 2 and a 3 ~ 0.2, this correction is of order 640 x e A £, 2 v R . 
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